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We examine the singularity resolution issue in quantum gravity by study- 
ing a new quantization of standard Friedmann-Robertson- Walker geometrody- 
namics. The quantization procedure is inspired by the loop quantum gravity 
programme, and is based on an alternative to the Schrodinger representation 
normally used in metric variable quantum cosmology. We show that in this 
representation for quantum geometrodynamics there exists a densely defined 
inverse scale factor operator, and that the Hamiltonian constraint acts as a 
difference operator on the basis states. We find that the cosmological singular- 
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ity is avoided in the quantum dynamics. We discuss these results with a view 



to identifying the criteria that constitute "singularity resolution" in quantum 
gravity. 
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m ', I. INTRODUCTION 

o 

q-i It is widely believed that a quantum theory of gravity will give insights into the question of what 
becomes of classical curvature singularities. This is based largely on intuition from uncertainty 



principle and fundamental length scale arguments in regions of large spacetime curvatures. What 
is required to address the problem quantitatively is quantization of model systems that contain 
classical metrics with curvature singularities. Such models are usually symmetry reductions of 
general relativity or other generally covariant metric theories. Within a model an obvious approach 
is to look at classical observables such as curvature scalars, and see if they can be represented as 
operators on a suitable Hilbert space. Their spectra and quantum dynamics may give an indication 
of what becomes of the classical singularity. 

This question has been studied using models derived from symmetry reductions of general rela- 
tivity since the late 1960 's 

aaaa. 

All of this work used the Arnowitt-Deser-Misner (ADM) 
(metric variable) Hamiltonian formulation of general relativity ("geometrodynamics") as the classi- 
cal starting point, and the Schrodinger representation as the quantum starting point for developing 
a quantum gravity model. The results obtained from various mini- and midi-superspace models were 
largely inconclusive. Some indicated singularity avoidance, others did not, but no insights emerged 
as general and definitive in the sense of transcending the model studied. 
After the development of the Ashtekar (connection variable) Hamiltonian formulation ("connection 
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dynamics"), many of the questions studied in the ADM formulation were revisited, including the 
general canonical quantum gravity program (for reviews see 0, 0] ) ■ The different classical variables 
led to the development of a non-Schrodinger representation program based on holonomy variables 
("loop quantum gravity"). Recently results from this programme were applied by Bojowald 0,0] 
to the old question of quantizing mini-super space models, with a view to studying what happens to 
classical curvature singularities upon quantization. This application has produced some interesting 
results: for Friedman-Robertson- Walker (FRW) mini-superspace models the Hamiltonian constraint 
acts like a difference operator on the space of states, and there is an upper bound on the spectrum 
of the inverse scale factor operator. Taken together, these results lead to the conclusion that the big 
bang singularity is resolved in the loop quantum gravity approach 0|. 

A number of questions may be asked at this stage concerning classical variables, quantization pro- 
cedures, and singularity resolution: What criteria constitute singularity avoidance? Is the singularity 
avoidance conclusion from the loop quantum gravity programme a result of both the classical start- 
ing point and the choice of representation? Would a non-Schrodinger representation quantization in 
the geometrodynamical ADM variables give the same results? 

Motivated by these questions, we study a new quantization of flat FRW cosmology (the model in 
which the loop quantum gravity results were first obtained 0,0]). Our classical starting point is the 
geometrodynamical Hamiltonian formulation, which we quantize via a non-Schrodinger representa- 
tion motivated by holonomy-like variables. We obtain results qualitatively similar to those obtained 
in loop quantum cosmology: the Hamiltonian constraint acts like a difference operator on the space 
of states, and the spectrum of the inverse scale factor has an upper bound. 

This paper is organised as follows: The next section introduces the classical theory and a set of 
basic variables for FRW geometrodynamics. Section III describes the quantization procedure and 
discusses the volume, inverse scale factor, and Hamiltonian constraint operators. The final section 
contains a discussion of the results with the singularity resolution question in mind. 

II. CLASSICAL THEORY 

The canonical Hamiltonian variables in geometrodynamics are {q a b,^ ab ), where the configuration 
variable q ab is the metric on a spatial 3-surface, and the conjugate momentum variable 7f ab is a 
function of the surface's extrinsic curvature. In terms of these variables, the vacuum Hamiltonian 
and spatial diffeomorphism constraints are 

H = (n ab n ab - ±tt 2 ) - = (1) 

C a = D b 7t b a = 0, (2) 

where q =detq ab , n = 7r ab q ab , and R is the Ricci scalar of q ab . 
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For the flat FRW model a reduced Hamiltonian theory may be obtained by the ansatz 

q ab = B(t)e ab , n ab = P(t)e ab (3) 

where e a b is the flat Euclidean metric diag(l,l,l). The density weight on 7f ab is obtained using this 
fiducial metric. Unlike in general relativity where Bit) is dimensionless since it is a metric variable, 
we will take it to have dimension length square (which means that the spatial coordinates are 
dimensionless). The conjugate momentum Pit) is then dimensionless in order that the symplectic 
form 

uj = -^—dBAdP (4) 
has dimensions of action (in c = 1 units). The phase space topology is that of a particle on the half 



line, R + x R, since Bit) > in this parametrization. The standard FRW scale factor is a(t) = y B(t). 
The Hamiltonian constraint is 

H = -\p 2 ^B (5) 

and the diffeomorphism constraint C a vanishes identically. It is interesting to note that this Hamil- 
tonian constraint is identical in form to that obtained in the connection variables Q. 

To facilitate quantization however, it is useful to make a canonical transformation that lifts the 
half line restriction. We proceed in two steps. The first is to extend the classical configuration space 
to include the singularity Bit) = 0. This step is essential for addressing the question of singularity 
resolution for without it there is "no classical singularity to avoid" in the quantum theory. (This for 
example is done - more or less unconsciously - in the familiar quantization of the hydrogen atom by 
the choice of Hilbert space L 2 (R^), even though the classical configuration space is only R + .) The 
second step is to reparametrize the configuration variable such that its domain is the real line. 

This is done by introducing new variables (x,p) by writing x 2 = B. The symplectic form becomes 

u = - - 2xdxAdP= o \ 1 dxAd{2xP). (6) 

Thus the new momentum is p = 2xP, for which {x, p) = 8ttGn- (Both coordinates x and p now have 
dimension length.) Note that in this parametrization, the point x = is included to give the full 
real line as the configuration space. This amounts to an extension of the original parametrization 
to include the degenerate metric with B(t) = 0. (This feature is also present in the connection-triad 
variables, where invertibility of the triad is not a requirement.) 
The Hamiltonian constraint ([!} as a function of (x,p) is 

H=--£\x\ = --£ (7) 
8 x 2 m 8 U 

Note that this constraint is now quite different in form than the one in the connection-triad variables, 
although the x variable may be regarded as a "triad." The reason for this is that the new momentum 
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p is a product of the metric and extrinsic curvature variables, and this is quite unlike the connection 



variable A made of the triad connection T and extrinsic curvature K as A = T + K. 



t is this 



IB- 



connection A that is essential for formulating the loop representation through holonomies 

We now introduce an algebra of classical observables and write quantities of physical interest 
as functions of these variables. Their form is motivated by the holonomy observables used for 
quantization in the loop quantum gravity programme. We use as the fundamental variables x and 

Uj(p) := exp(i7p/L), (8) 

where 7 is a real parameter and L fixes the (arbitrary) unit of length. The parameter 7 is necessary 
in order to separate momentum points in phase space (eg. fixing 7/L = 1 say, gives the same value 
of U for p and p + 2mr). This variable may be viewed as a momentum analog of the holonomy 
variable of loop quantum gravity. 
The pair (x, U^ip)) has the Poisson bracket algebra 

{x,U 1 (p)} = (8nG N ) l -l U 7 (p). (9) 

This is the basic algebra which will be taken over to the quantum theory. 
Other quantities of interest are the volume, which up to a multiplicative constant is 

V(t) = B{tf' 2 = |x| 3 , (10) 

and inverse powers of the scale factor a(t). The latter has proven useful in determining whether 
there is "singularity avoidance" in a quantum theory |2|. It is possible to write these observables, 
and the Hamiltonian constraint using our choice of fundamental variables (x, U^{p)). 
The standard FRW scale factor is 



a (t) = y/B(t) = \x\. (11) 

Using the method introduced in [l(]], inverse powers of the scale factor may be defined classically 
either via the bracket 

U- x {V n ,U^} = ^{M 3 ",^} = i{MG N ) J 3n sgn^lxl 3 "" 1 , (12) 

or more simply by inverse powers of the volume observable. However the latter definition cannot be 
carried over to the quantum theory since, as we see will below, the volume operator turns out to 
have zero as an eigenvalue, so negative powers of it are not densely defined. Therefore the somewhat 
indirect definition (J12|) . with n > will be useful in studying the approach to the singularity in the 
quantum theory. The requirement that the power of x on the right hand side be negative means 
that n < 1/3. Thus we need 0<n<l/3in order to obtain a sequence of inverse powers of the scale 
factor in terms of the basic variables. The choice n = 1/6 gives 

sgn(x) 2Li 



(SvrGVh 



U^iV 1 /",^}, (13) 
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and n — 1/4 gives 

1 



(f/- 1 ^ 174 ,^}) 4 (14) 



M V3(87tGjv)7L / 

From this it is clear that powers of Poisson bracket in Eqn. (fT2|) may be used as a starting point 
for defining a large class of operators for inverse powers of the scale factor. It is of interest to see 
whether these all lead to qualitatively similar behavior concerning the quantum nature of the big 
bang singularity. 

The Hamiltonian constraint can be written as a function of the basic variables by using the relation 
f|13|) for the inverse scale factor as follows: 

^irw' 2 ^'''^ (15) 

Note here that there is the alternative choice of using (fPfj) to define the Hamiltonian constraint. 
However, this leads to a more complicated form due to the larger number of U factors. 

It is of course also possible to write the Hamiltonian constraint in other classically equivalent ways. 
One alternative is to substitute \x\/x 2 rather than directly using the inverse scale factor l/\x\ from 
f!13|) or (jl4j) . These choices will clearly lead to inequivalent operators in the quantum theory since 
the number of factors of U are different. As for the inverse scale factor, it is generally useful to also 
study various choices for the Hamiltonian constraint in order to identify the main features that are 
to be associated with "singularity resolution" in the quantum theory. 

In the following, we focus primarily on the simplest ordering of the Hamiltonian constraint, and 
the square of the definition (fTH|) for the inverse scale factor, since both of these contain the smallest 
number of U factors. However we will briefly comment on other choices. 



III. QUANTUM THEORY 



To construct the quantum theory for the classical system described above, we will proceed in 
analogy to the procedure used in loop quantum gravity. The first step is to choose an algebra of 
classical functions that is represented as quantum configuration operators. We take here the algebra 
generated by the functions 

W{\) = e lXx/L , (16) 
where A6K. It consists of all functions of the form 

n 

f(x)=Y / c,e iX ^ L , (17) 

with Cj G C and their limits with respect to the supremum norm. This algebra is known as the 
algebra of almost periodic functions over R and we denote it by AP(R). 
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It is well-known that AP(R) is naturally isomorphic to C(W.Bohr), the algebra of continuous func- 
tions on the so-called Bohr-compactification of K ^lj. As the name suggests, ^Bohr is a compact 
group which can be obtained as the dual group of Rd,iscr, the real line endowed with the discrete 
topology. This suggests that taking L 2 (WBohr,dfio), where /xo is Haar measure on Rsohr, as the 
Hilbert space for our theory is a viable option. This is the decisive point where we depart from the 
traditional approach in geometrodynamics, where the Hilbert space is the conventional Schrddinger 
space L 2 (R, dx). Once we adopt this new choice, basis states in our Hilbert space are given by 

| A) = \e iXx / L ), A e M, (18) 

with the inner product 

</i | A> = <W- (19) 

This representation has been discussed in some mathematical detail in 0], and also in where it is 
applied to the quantization of a particle. Notice the difference from the standard quantum mechanics 
of a particle on the real line, where the right hand side is given instead by a delta function 5(fi — A). 
This feature is traceable to the fact that the configuration space is the real line with the discrete 
topology, which in turn stems from the choice of the algebra of functions. 
The action of the configuration operators VF(A) is defined by 

W(\)\/j) = e iA£/L |/i) = e^|/i). (20) 

It is straightforward to verify that these operators are weakly continuous in A, which procures the 
existence of a self-adjoint operator x, acting on basis states according to 

x\n) = Ln\n). (21) 

The next step is to construct the operators corresponding to the classical momentum functions 
U y = e* 7P//L . Their action on the basis states is fixed by the definition of the x operator and the 
requirement that the commutator between x and Uj reflects the corresponding Poisson bracket Q 
between x and U T With the definition 

£» = |/i- 7 >, (22) 

the commutator is 

[x, U 7 ) = -jLU^ (23) 

Making now the standard commutator- Poisson bracket correspondence [, ] < > ih{, }, gives using 

(0) the relation 

-7.L = ih (8vrGjv)^, (24) 
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which fixes the length L to L = y&Klp. This shows explicitly how the eigenvalues of x arise in 
Planck units upon quantization. 

Obviously, ?7 7 is unitary, however, it is not weakly continuous with respect to 7. As a conse- 
quence, there is no momentum operator in this representation, in stark contrast to the Schrodinger 
quantization. 

With the basic quantum operators now at our disposal, we are in a position to construct the inverse 
scale factor operator and investigate its spectrum. 

A. Volume and Inverse Scale Factor 

The operator for the volume V is provided directly by the operator x defined in ()21j1 . We have 

V\n) = (v^/p) 3 !^! 3 ^). (25) 

The operators corresponding to U and V can be used to obtain an operator for the inverse scale 
factor. One way to do to this is to use the square of the expression in Eqn. (|THj) with 7 = 1. The 
resulting operator is 

The key question is whether this operator is unbounded as in standard quantum cosmology, where 
its eigenvalues diverge when approaching the quantum state corresponding to a = 0, or whether 
it is bounded, indicating a (kinematical) resolution of the classical singularity. To decide this we 
calculate its eigenvalue on a basis state 




(27) 

This result reveals some important properties of the eigenvalues. First, they are always positive or 
at most zero, as should be the case. Second and more importantly, the spectrum is clearly bounded 
from above. For \fi\ — > 00 the eigenvalues approach 0, as would be expected from the behavior 
of l/\x\ for large Moreover, the eigenvalue of the state |/i = 0) corresponding to the classical 
singularity (a|0) = |x||0) = 



is y2/7r/p, and this is the largest possible eigenvalue. (This is notably 
different from the results in [9], where the eigenvalue of the inverse scale operator for the state \fi = 0) 
is 0, and the maximal eigenvalue is obtained instead for the state \/i — 1). Although there are no 
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principal reasons why this could not happen in the quantum regime, it seems somewhat unnatural 
from the classical point of view. It should be pointed out however that this result is obtained in our 
formalism for a different choice of operator ordering.) 

In summary, this new quantization of the inverse scale factor in geometrodynamics mimics the 
expected classical behavior for large values a(t), and departs significantly from the divergence in the 
standard quantization near the classical singularity a(t) = 0. In this sense, the quantization resolves 
the singularity. This "resolution" however is so far only kinematical, since we have not investigated 
the quantum dynamics. It is conceivable that the quantum dynamics breaks down at the state |0), 
in which case it would be hard to claim a satisfactory resolution of the singularity. As the dynamics 
is encoded in the Hamiltonian constraint, we now turn our attention to its operator realization. 



B. Hamiltonian constraint 



As discussed already in the classical section, the Hamiltonian can be written in many different, 
classically equivalent forms. The one we will focus on in this section is 

(28) 

as this is in some sense the simplest one, and the spectrum of the inverse scale operator is already 
known. As p does not exist as an operator in our quantum representation, we have to choose an 
alternative way to represent p 2 as an operator. One way to do this is motivated by the classical 
expression 

pWlim 1 (2 - U, - U?) • (29) 

A physical interpretation of this expression is obtained by setting 7 = lp/L p h ys where L p h ys is the 
characteristic size of the system under consideration, and lp is a fundamental length scale. (Note 
that a Hamiltonian naturally introduces a scale L p h ys for a physical system.) The limit then suggests 
that the "point" form of the momentum is recoverable in the case L p h ys » If- 

For quantum cosmology these considerations mean lp = lp and 7 = lp/L p h ys , and lead to a 
Hamiltonian constraint operator 

tf 7 = ?< (t> 7 + £>-!_ 2 ) J 

= ^fr + u; 1 -2)(v- 1 \ywv])\ (30) 

where a specific operator ordering has been chosen. The action of on a basis state is given by 



HM = ^/p(H 1/2 -||i-ir /2 ) 2 (|^ + 7> + l^-7>-2||i» 



r 
I 8 

7 2 



IpV(fx) (|/i + 7 > + |/i-7>-2|^» (31) 
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On the eigenstate |0) of volume with zero eigenvalue, which is the classical singularity, we have 

H,\0) = ^ l P (h) + I - 7) " 2|0)), (32) 



R l°> = )/4 |0> (33) 

These equations represent the effects of quantization on the classical singularity. In order to probe the 
dynamical part further we must solve the quantum constraint equation that encodes time evolution. 

As is well known in the theory of constrained systems, normalizable solutions of the quantum 
constraints do not lie in the kinematical Hilbert space 7i, but rather in a larger space C*. This space 
can be obtained as the dual space of the dense subspace C of Ti, which is spanned by all elements of 
the form 



5>aoi^>- (34) 
i=i 

A general element of C* can thus be written as 

(VI = £VWM- (35) 

Notice that, while the sum is continuous as it runs over every real number, its action on an element 
of C is well defined by construction. The constraint equation - symbolically written as 

H\i>) = 0, (36) 

is now interpreted as an equation in the dual space, 

(V>|# f = 0. (37) 

Using the form of a general element of the dual space (|35j) and the (dualized) action of the (dual) 
Hamiltonian on (dual) basis elements, we can derive a relation for the coefficients ip(fi): 

V(ji + 7 M/i + 7) - 2V(/iM/i) + V(/i - 7 M/i - 7) = 0. (38) 

What is the meaning of this equation and in what sense does it encode the quantum dynamics? 
First of all, it determines the coefficients for those dual states that are physical. As in the classical 
theory solutions to the constraint equation represent classical spacetimes, these physical dual states 
can be interpreted as representing " quantum spacetimes" . 

The difference equation (|38|) gives physical states are linear combinations of a countable number 
of components of the form 

^(ji + n^lfi + wy), (39) 

where 7 is fixed at the Planck scale (7 = lp/L p h ys ~ 1) and jiGZ. As each component corresponds 
to a different eigenvalue for the volume and scale factor, it can be interpreted as the quantum 
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state representing the universe at the "time" \i + 7. A solution of the Hamiltonian constraint 
therefore represents a linear combination of FRW universes specified at certain discrete volumes, or 
equivalently, at discrete times. It is in this sense that time evolution is "discrete with fundamental 
time step" 7. It is also clear that this "discrete evolution" does not represent the state of a single 
universe at different discrete times, since the term "single universe" has no meaning here. Rather 
a "discrete solution" of the Hamiltonian constraint, (ie. one satisfying (fHKJ) ). gives the amplitudes 
that the physical universe is in one or other of the discretely separated components of the physical 
state. 

The state < 0| corresponding to the classical singularity is contained in only one specific "quantum 
spacetime," (ie. solution of the Hamiltonian constraint). Furthermore, in that one case we can see 
that the system evolves right through the singularity without encountering any problems, since the 
component ip(0) can be computed in terms of the components V'(t) an d 7)- In all other physical 
states, the state < 0| does not occur, and so in a sense one can say that the discrete evolution 
"jumps" over the singularity if the state contains components with both positive and negative values 
of \x. In such cases there is an instance of smallest but finite volume. 

From these observations one can perhaps conclude that dynamically the singularity has been 
resolved. A dynamical non-resolution of the singularity might have occurred had it turned out that 
the difference equations have no solutions if they contain the ip(0) component, or if they contain 
components with both positive and negative /1 in the sum (J35|) . 

Finally, it is interesting to note that for our representation of H, the state at the classical singu- 
larity "^(O) can be denned in contrast to the results in & However, had we chosen to write 
the classical Hamiltonian using eqn. (jl4j) instead of eqn. ()13|). which amounts to using double the 
number of U operators, we would have ended up with the same result: ^(0) cannot be determined 
from the difference equation, but a solution is still possible as it turns out that VKt) is then given 
in terms of 7). This shows the significant differences that can arise due to quantization ambigu- 
ities. Ultimately, only physical predictions and comparison with know facts or (as yet hypothetical) 
experiments can determine the "right" choice. 

IV. CONCLUSIONS AND DISCUSSION 

Our main result is that there is an alternative to the Schrodinger quantization of the FRW cos- 
mology in the standard ADM geometrodynamical variables. This quantization leads to conclusions 
qualitatively similar to those obtained in loop quantum cosmology starting from the connection- 
triad variables: (i) the Hamiltonian constraint acts like a difference operator, and (ii) the inverse 
scale factor can be represented as a densely defined operator. Thus it is the representation space 
and the realizations of the basic observables rather than the nature of the classical variables that are 
responsible for the similar conclusions for this model. 

To what extent is the quantization we have presented different from the one employed in loop 
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quantum cosmology? The differences at the classical level are clear: the phase space variables 
(x,p) are not the standard mini-superspace variables that arise via standard reduction from the 
connection-triad canonical variables, as comparison with 0, 0| shows. The key difference at the 
quantum level is that U 1 is not the holonomy operator associated with the Ashtekar-Sen connection 
for the FRW model. Rather, the LL we use is a standard translation generator whose realization on 
the Hilbert space L 2 (M.Bohr,d[io) is applicable to any classical theory, as has been discussed in [12]. 
Thus interpreting our quantization as a "loop representation" would mean that one is generalizing 
this term to include all quantizations on the Hilbert space L 2 (M.Bohr,dfi ). 

It is clear that the alternative representation based on the Bohr compactification is applicable to 
other mini-superspace models, since in all such models the phase space variables are functions of only 
a time coordinate. It is also clear that this applicability is independent of whether the classical phase 
space variables are metric-extrinsic curvature or connection-triad. The main difference between the 
variables arises in the form and action of the Hamiltonian constraint. 

In the flat FRW case we have discussed, the Ricci scalar term in the Hamiltonian constraint 
vanishes. Thus the action of the constraint as a difference operator is due only to the kinetic term. 
In other mini-superspace models the Ricci scalar term, which is a purely configuration variable, 
will have non-trivial action on the basis states. However in the Bohr representation, this action 
is multiplicative. Thus it appears that in other mini-superspace models the "difference operator" 
feature of the Hamiltonian constraint will survive. Similarly it appears that an inverse scale factor 
operator is definable using volume and U operators, and that it is likely to have a spectrum bounded 
above. For models where the phase space is more than two-dimensional, the new representation can 
clearly be used for each pair of phase space variables. Extension beyond mini-superspace (quantum 
mechanics) to midi-superspace (quantum field theory) models, such as the Gowdy cosmology would 
be of much interest 14 1. 

"Singularity resolution" appears to consist of two main features, one kinematical and the other 
dynamical. The kinematical feature is the spectrum of the operator associated with a curvature scalar 
(or other relevant classical observable) that diverges at a curvature singularity. If the spectrum is 
bounded, the singularity may be considered kinematically resolved. It is important to identify the 
largest eigenvalue and corresponding eigenstate of such an operator, since this is the "closest" the 
quantum theory can get to the singularity. The dynamical feature of singularity resolution concerns 
the action of the Hamiltonian constraint on the state of largest curvature: this could lead either to 
no solution of the constraint for zero or negative values of /i, or to a well defined "evolution" through 
zero to negative values of [i. The former may be taken as an indication of the breakdown of quantum 
evolution, and hence a dynamical non-resolution of the singularity, regardless of the boundedness of 
the curvature operator. 

An alternative viewpoint is that the kinematical vs. dynamical views are artificial in that the 
question of singularity resolution is relevant only for the physical state space with a well defined 
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physical inner product. The question then becomes whether there are any physical states for which 
the curvature operator spectrum is unbounded. This appears more compelling, but it has not been 
addressed here, or in the context of loop quantum cosmology. 
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